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Abstract 

We classify all the six derivative Lagrangians of gravity, whose traced field equations are of 
second or third order, in arbitrary dimensions. In the former case, the Lagrangian in dimensions 
greater than six, reduces to an arbitrary linear combination of the six dimensional Euler density and 
the two linearly independent cubic Weyl invariants. In five dimensions, besides the independent 
cubic Weyl invariant, we obtain an interesting cubic combination, whose field equations for static 
spherically symmetric spacetimes are of second order. In the later case, in arbitrary dimensions 
we obtain two combinations, which in dimension three, are equivalent to the complete contraction 
of two Cotton tensors. Moreover, we also recover all the conformal anomalies in six dimensions. 
Finally, we present the general static, spherically symmetric solution for some of these Lagrangians. 
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1 Introduction 

Einstein's General Relativity is not only the most successful classical theory of gravity in four dimen- 
sions, it is also the simplest theory possessing some nice properties. The two of the most important 
characteristics being general covariance and second order field equations. In fact, it was shown by 
Lovelock that, in four dimensions, General Relativity is the unique generally covariant theory of grav- 
ity (upto an addition of a cosmological constant) which gives second order equations of motion [T]. 
However, in higher dimensions, there exists higher curvature theories, namely the Lovelock theories, 
which also gives second order field equations. These theories are generically characterized by higher 
curvature invariants in the action. At each order k, the combination of higher curvature invariants 
is unique, the integral of which on a compact manifold of dimension 2k gives the Euler characteristic 
of the manifold. There are also other interesting characteristics of the Lovelock class of theories. 
Primarily, exact analytic black hole solutions are known to exist [2j [3] (see also [U [5] and references 
therein). Black holes are widely believed to exist in nature as a final state of a gravitational collapse. 
They are also the simplest objects to study in any gravitational theory. Therefore, exact black hole 
solutions are of significant importance. More recently, in the context of AdS/CFT correspondence, 
exact asymptotically AdS black hole solutions in a gravity theory have been proven to be useful in 
studying the holographic properties of a finite temperature conformal field theory on the boundary [6]. 
Secondly, the Lovelock theories also admit Birkhoff's theorem, which states that any solution which 
has spherical, planar or hyperbolic symmetry must be locally isometric to the corresponding static 
black hole solution [7]. Generically, the admittance of Birkhoff's theorem suggests the lack of spin-0 
mode excitations in the linearized field equations. 

There are also other theories of gravity which admit exact analytic black hole solutions and fur- 
ther admit Birkhoff's theorem. These theories, being outside the Lovelock class, are generically higher 
derivative theories and consequently possess ghost degrees of freedom. One well known higher deriva- 
tive theory is the conformal theory of gravity in four dimensions, which is obtained from an action 
quadratic in the conformal Weyl tensor. The action is thus invariant under Weyl rescalings and the 
field equations are traceless. Birkhoff's theorem in conformal gravity states that modulo a conformal 
factor the most general spherically symmetric solution is static [8]. Considering the same action in 
dimensions other than four, one loses the property of invariance under Weyl rescalings. Even then, the 
theory admits exact analytic black hole solutions with spherical, planar or hyperbolic symmetry, and 
further admits Birkhoff's theorem [9]. Note that, in arbitrary dimensions, though the field equations 
are of fourth order, the trace of the field equations are of order two. This can easily be seen as follows. 
Varying the action gives 



Now, consider infinitesimal Weyl-rescalings of the metric 5g^ u = \g^ v '. Under such variations, the 
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Lagrangian will vary as A(2 — D/2)^/—g£, which implies = (^— )£. This shows that the trace of 
the field equations, being proportional to the Lagrangian density, must be of second order. 

Recently, a theory of massive gravity was constructed in three dimensions [lOj . where the La- 
grangian is a particular combination of quadratic curvature invariants given by 

K = R ab R ah - Z -R 2 . (2) 

Again, the theory admits exact analytic black hole solutions [llj. The Lagrangian of this theory has 
a unique property in three dimensions, that the field equations have a second order trace. So, it is 
natural to wonder if there could be other theories of gravity which, although non-realistic, have some 
special properties which allows one to obtain exact analytic black hole solutions and can serve as toy 
models of gravitational theories. Specifically, it might be useful to classify higher-derivative theories 
of gravity whose traced field equations have a reduced order. 

In this paper we construct the most general Lagrangian which is a linear combination of scalars of 
the form 

/»'.... /»'.... /.'.... and V.22....V.J?.... , (3) 

which are characterized by the number of derivatives of the metric (hereafter, the degree of differen- 
tiation) n = 6, such that the trace of the field equations is of order three or less. 

We will show that when the trace is resticted to be of order two, then, in dimensions six or 
higher, there are only three linearly independent possible invariants which have a second order trace. 
They are the six-dimensional Euler density, and the two linearly independent scalars constructed by 
contracting all the indices of three conformal Weyl tensors. However, in five dimensions, we obtain a 
peculiar independent invariant which can be thought of as a cubic generalization of ([2]) and has been 
studied separately in [12] and [131 HI] • We also obtain the general static spherically symmetric solution 
for some of these theories. Based on our analysis for six-derivative theories, we present a conjecture 
classifying all the scalars of arbitrary order, which give second order traced field equations in various 
dimensions. 

When the trace is restricted to be of order three, in arbitrary dimensions, we obtain two additional 
scalars. These two scalars are not independent in dimensions three and six. In six dimensions they 
reduce to one of the conformal anomalies. 

One future direction of study is to see which of these theories admit a Birkhoff 's theorem. 

In Section [21 for completeness we review the n = 4 case. The case n = 6 is analyzed in Section [3l 
In Section H] we focus our attention on obtaining the general static, spherically symmetric solution for 
some of the theories defined in Section [3] in arbitrary dimensions. 

2 Quadratic Combinations, n = 4 

In this section, we review how to obtain the most general quadratic Lagrangian, having second order 
traced field equations [HI [16] . 
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The most general quadratic combination of curvature invariants in arbitrary dimension is given 



It 
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C Q := a R abcd R abcd + b R ab R ab + c R 2 , (4) 



where a, b and c are arbitrary constants. The trace of the field equations coming from this Lagrangian 
are 

G^ a = ^4a + ^b + 2 (D - 1) cj V a V a R - AaV a V b R ab - ^~j^£q- (5) 

Imposing the trace G^ a to be of second order, implies that the coefficients a, b and c must be chosen 
such that the first two terms in ([5]) vanish, i.e., 

D 



V„V b 



AaR ab + g ab [Aa + ^b + 2(D - l)c ) R 



0. (6) 



The above equation is satisfied only if the term inside the bracket is proportional to the Einstein 
tensor, which is the most general divergenceless, symmetric, rank two tensor linear in the curvaturj§. 
Consequently, the coefficients in ([6]) must fulfill 

-4a = 7 and 4a + ^ b + 2(D - l)c = -~ . (7) 

Since there are four variables (a, 6, c, 7) and two equations, there is a bi-parametric family of solutions 
given by 

C Q = -lR abcd R abcd + b R ab R ab + ^^ R 2 ■ (8) 
One can further factor out the four dimensional Euler density to write ([8]) in the form 

C Q = aN± + PS A , (9) 
where a = b — 7 and /3 = — j are arbitrary constants, A/4 is defined by 

\r Aj? ab i? ^ p 2 — ^ 1 ^ ~ ^ 1 x a i b i a 2b 2 ( r< c\d\r< c 2 rf 2 p p c 2 d 2 \ / ln \ 

A/ 4 .- 4rt K ab - Jp—[) K ~ 2? \ D^3 J c i rf i c 2^2 [^aih C a 2 fe 2 ~ ^axW K a 2 b 2 ) ' l iU J 

and £4 := i? 2 — AR ab R ab + R a bcdR abcd is the Gauss-Bonnet combination which corresponds to the four 
dimensional Euler density. Thus, we have shown that for combinations quadratic in the curvature, the 
most general Lagrangian which has second order traced field equation, can be expressed as a linear 
combination of the four dimensional Euler density £4 and A/4 defined in Eq. (jlOp . 
In dimensions higher than three, one can further use the following relation 

C abcd C abcd = £ 4 + A/4 , (11) 



Note that the only non-quadratic term with degree of differentiation 4 is OR, which is boundary term. 
2 Since a divergenceless vector J a cannot be constructed locally out of the curvature, the equation V a J a = 0, with 
J a := V 6 [~4aR ab + g ab (4a + §-b + 2(D - f)c) 7?] , does not have any non-trivial solution. 
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where C abc d is the Weyl tensor and C C abc d is the four-dimensional conformal Weyl invariant. This 
means that for dimensions higher than three, Cq in ([5]) can be equivalently expressed as a linear 
combination of the four-dimensional Weyl invariant and the Euler density. 



3 Lagrangians with n = 6 

In this section, we generalize the previous discussion for Lagrangians of degree of differentiation six 
in arbitrary dimensions D. We start by considering a generic combination of the twelve linearly 
independent |18j . curvature invariants of degree six 

12 

£ = J>%, (12) 

i=i 

where ^4*'s are arbitrary coefficients and Li's are given by 

r riabcd p p e / r jyab pee pdf j r>abcd p pe j p ryabcd p j r>abcd p p 

L\ — it H cde fH ab , L 2 — K cd H b jK ae , L 3 — K H cdbe K a , — titi ti ahcd , L 5 — H ti ac ti hd , 

Lq = R ab Ri )C R c a , L-j = RR ab R a }j, Lg = R 3 , Lg = V a RV a R, Liq = V a Rb c V a R bc L\\ = V p R a b cd \7 p R abcd , 

L12 = V a R bc S/ b R ac . (13) 

Note however that when we neglect a total derviative the twelve terms are not linearly independent, 
as one can write two of the invariants in terms of the other ten in the following waj§ 

L u = 2Li - 4L 2 + 2L 3 - 4L 5 + 4L 6 - L 9 + 4L 10 + V a V c [2R abde R c bde - 8R abcd R bd + 8RR ac 
- 12R ab R c b + 2g ac (2R bd R bd - R 2 )} 

Ln = L 5 -L 6 + ^L 9 + V a V c [^ a 6J R bc - RR ac + \g ac R 2 } (14) 

Using the above relations, one can rewrite the Lagrangian as a linear combination of only ten curvature 
invariants with coefficients A 1 {i = 1, • • • , 10). Extremizing the action constructed with the Lagrangian 
(fT2l) with respect to the metric gives the field equations 

10 

^=£^[^ = 0, (15) 

i=l 

where G/J b are defined respectively in equations (|45H56p in the Appendix |Bl Now, requiring the trace 
Qa^ a to be of order three, is the same as imposing trace to be proportional to C (See Appendix |A|) . 
This gives us a set of eight equations (|58p for the twelve variables (^4*) and a parameter u (analoguous 
to 7 in the quadratic case (O). Solving these equations for arbitrary dimensions, we obtain a five- 
parameter family of solutions. The details of the equations and its solution are given in the Appendix 
IB] This implies that in D > 6, there are five linearly independent curvature invariants (of degree 



We would like to thank Nicolas Boulanger for pointing this out to us. 
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of differentiation six) which gives rise to third (or lower) order traced field equations. They are as 
follows. Firstly, the six-dimensional Euler density given by 



£q := 2L X + 8L 2 + 24L 3 + 3L 4 + 24L 5 + 16L 6 - 12L 7 + L 



8 • 



(16) 



obviously gives second-order traced field equations. Secondly, there are two independent algebraic in- 
variants constructed out of three Weyl tensors, namely W\ = C ab cd C cd e jC e ^ ab and W 2 = C a b c dC ebc ^C a e j- d , 
which also gives second order traced filed equations. These two Weyl invariants are given in terms of 
the L,-'s as 



W 1 = L 1 + 



12 



+ 



D - 2 
16(D - 1) 



L 3 + 



6 



(D - 2)3 



(D-l)(D-2 
2A(2D - 3) 
(£>-!)(£> -2)3 



La + 



24 



L 7 + 



{D-2) 2 ' 

8(2D - 3) 
(D- l) 2 (D-2) 3 



I* 



(17) 



and 



Wo 



1 3 
- A Li + L 2 + — 



-L 3 + 



2(D - l)(D - 2) 



La + 



3D 



(D - 2) 2 



2{3D - 4) 
+ (D - 2)3 Le 



3(D 2 + D-4) 
(D-l)(D-2f 



L 7 + 



(D 2 + D- 4) 
(L»-l) 2 (D-2) 3 ' 



(18) 



In addition, there are two other curvature invariants £ and listed in equations (|6ip and (|62p 
respectively which give third order traced field equations. However, in dimensions D < 6, the above 
curvature invariants are not all linearly independent. For example, in D = 3 and 6, the invariants 
X and are proportional to each other modulo a total derivative. It is interesting to note that, in 
six dimensions, requiring the traced field equations to be of third order (or less), we recover all the 
four (1 type-A and 3 type-B) non-trivial conformal anomalies [I2]-[5T]. In Table I below, we list all 
the curvature invariants in dimensions greater than or equal to 3, which lead to third (or less) order 
traced field equations. 



G a a 


D = 3 


D = 4 


D = 5 


D = 6 


D > 6 


d' 2 g 


$ 


Wi ~ w 2 


Wi ~ w 2 , Mq 


W h W 2 


W U W 2 , £ 6 




ri riabc 


s, e 


s, e 


s~ e 


s, e 



Table 1: Here C a b c denotes the Cotton tensor. Note that in dimension five a new combination 
Af 6 := -24L 3 - %L 4 - 40L 5 - 3|o Lg + ot^ _ ^ Lg appears . 

Returning to the set of invariants that gives second order traced field equations we find that in 
dimensions they are spanned by the basis set W±, W 2 } up to a total derivative. However, in 

D = 5 this is not the case. In particular, there exists a "special" linearly independent invariant which 
generalizes A/4 to the cubic case. This is realized by noting that the following relation is analogous to 
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equation (fTTj) 

4Wi + 8W 2 = £ 6 + (^t-t; ) ■A/e (19) 



where 



3(2) + 2) 24L> 16D(D-1) 12(D 3 - 2D 2 + 6D - 8) 
:= " 24Ls " (D - 1) L4 " ^2 Ls " (2? -2)2 L6 + (D — 2) 2 (D — 1) L? 
(D 4 - 3D 3 + 10D 2 + 4£> - 24) 
(2? -2)2(2? -1)2 Ls (20) 

is the cubic counterpart of A/4. Let us rewrite equation (fT9|) in the form 

M 6 := ^4(4^ + 8W 2 - £ 6 ) (21) 
V — 5 

— J_ 1 _ ~ ^ 1 X«l ^112^2 013^3 / ci^ir 1 c 2d2^< c 3 d 3 _ p cidijy c 2 d 2 p c 3 <2 3 \ /r> 9 \ 

~~ 2 3 \D - 5/ c i rf i c 2^2C 3 <i3 ^ U a!fe! u a 2 6 2 U a 3 6 3 ^0x61 ""aafea ^363 ; ^> 

The term inside the parenthesis on the right hand side vanishes identically in dimensions lower than 
five, since for D < 5 

£ 6 = 4T^i + 8W 2 = 0. (23) 

However, in D = 5 (and greater than 5) this gives a non- vanishing invariant as can be seen by 
expressing W\, W2 and £q in terms of {Li,-- - , I/g}, thereby obtaining the expression ([2"Dj) . This 
imples that in D = 5, the basis is {W\ (~ W2) jA/^} up to a total derivative. Similar results have 
been found for quartic invariants (see Appendix B of reference [12]) i.e., in dimensions D 7^ 7, any 
invariant giving second order traced field equations can be expressed as a linear combination of the 
eight-dimensional Euler density, and all the linearly independent Weyl invariants, however, in D = 7 
there is an additional "special" invariant which completes the basis. Based on these results, we present 
the following conjecture 

Conjecture: (i) In dimensions D ^ 2k — 1, any curvature invariant of order A;@ ; which gives sec- 
ond (or less) order traced field equations can be expressed as a linear combination of the 2k -dimensional 
Euler density, the Weyl invariants and a total derivative. 

(ii) In dimensions D = 2k — 1, any curvature invariant of order k, which gives second order traced 
field equations can be expressed as a linear combination of the Weyl invariants a total derivative 
and a "special" invariant which can be obtained by evaluating 

\f .— J_ ( D — 2 \ r- ai bi---a k b k ( f-i cidi f~i c k d k _ p cidi p c k d k \ / a\ 
^2k ■- 2k y D _2k + l) °cidi-c k d k y^aib! ' ' ' °a fe fe fe ^oifti " " " U a k b k J 

in D = 2k- 1. 



4 By a curvature invariant of order k, we mean a scalar constructed out of k curvature tensors without any derivatives 
acting on them. 

5 Note that the number of linearly independent Weyl invariants of order k in dimensions D = 2k — 1 is one less than 
that in dimensions D > 2k, due to the identity C, . . . C ? k \ bk = 0. 
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We now show that Mzk evaluated in D = 2k — 1 indeed gives second order traced field equations. 
First, consider the following invariant of order k 

1 raibi---a,kbk ( (~t c^d\ p c k d k _ p cidi p c k d k \ /«r\ 
2 k°c 1 d 1 -c k d k \ ( - / a 1 b 1 "'^a k b k ^ai&j " ' JX a k b k ) 

Obviously, the above invariant vanishes in dimensions lower than 2k. However, if one expands the 
Weyl tensor in terms of the Riemann tensor, then it can be factorized by (D — 2k + 1). This can be 
seen as follows. Consider the basis set of k-th order Riemann invariants in arbitrary dimensions. In 
D = 2k — 1, not all elements of this set are linearly independent. In fact, the basis set contains one less 
invariant than in D > 2k. This is because of the vanishing of the fc-th order Lovelock density. Now, 
after the expanding in terms of the Riemann tensors, the term (|25[) will not contain any (Riemann) k . 
So, this invariant cannot vanish identically in D = 2k — 1 unless it is factorized by (D — 2k + l)j§ 
Further expanding all the Weyl tensors, one can convince one self that the dimensional dependence of 
the coefficient of the term with k— 1 Riemann tensors and one Ricci tensor must be (D—2k+l)/(D—2). 
We can now divide this factor out to get a non- vanishing invariant in D = 2k — 1. Thus, we write the 
A;th order generalization of A/4 by evaluating 

J_ [ D — 2 \ r- ai bi---a k b k (p c\d\ f~i c k d k _ p cidi p c k d k \ /r^x 

2 k \D -2k+l) cldl "' Ckdk \ aibl " a kb k *V&i '"^a k b k ) 

in D = 2k — 1. Note that, by construction, the trace of the field equation arising from the above 
invariant is of second order in all dimensions. 



4 Exact Solutions 

In this section, we present exact, static solutions for the theories defined previously. For simplicity 
we will first focus on the theories having fourth order field equations, defined by an arbitrary linear 
combination of the invariants W\ and W2 , defined respectively in equations (fTTJ) and (fTBj) . The theory 
defined by the combination A/6 in the Table I, has further interesting properties in five dimensions, 
which we discuss in detail in reference [12] . Finally we comment on the new three dimensional theory 
shown in Table I, which possesses third order traced field equations. 
The class of metrics considered is: 

jp2 

dsl = -F (R) dt 2 + ^- + R 2 dY? D _ 2>1 , (27) 

where eE£>_2 i7 is the line element of a (D — 2)-dimensional compact, orientable Euclidean manifold 
of constant curvature 7. For 7 = 1 the manifold X_d_2 is locally equivalent to the sphere S D ~ 2 , while 
for 7 = it reduces to a locally flat manifold. Finally for 7 = — 1 the geometry of S^>_2 is given by 
the quotient Hd-2/F, where T is a freely acting, discrete subgroup of O (D — 2, 1). 

6 This argument cannot be extended to dimensions 2k — 2 since one obtains another identity involving the Riemann 
invariants which is obtained be contracting the Ricci tensor with the (k — l)-th order Lovelock equation. 
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After a coordinate transformation and a redefinition of the arbitrary functions, the line element 
d2Z|) takes the form 

r A** . _ i 

(28) 



ds^ = N (r) 



As shown below, this gauge choice is much more convenient for our purposes. 
4.1 C 3 theories 

Here we will consider Lagrangians of the form 

£ = aWx + fiW 2 . (29) 

It has been proved in [22] that, for the metric (|28p . the two invariants W\ and W% defined in (|17p and 
(|18|) respectively, are proportional. Consequently, for a particular choice of at//3, both the Lagrangian 
C and the field equations vanish identically. In such a situation any metric within the family (|28p is 
a solution of the system. Hereafter we assume that a and (3 are generic. 

Since, in six dimensions the gravity theories defined by combinations of W\ and Wi are invariant 
under Weyl rescalings, let us concentrate on this case first. 

D = 6 : Using Weyl rescalings, one can gauge away the function N (r) in (|28|) . Then, one finds 
that the solution reduces to 

dsl = - (ar 2 + br + K- dt 2 + dl * - ^ + r 2 dY?^ , (30) 



/ ar A + or + .« 

where a, 6, AT, c and e are constants, which are related by 

A = 7 and c(b- 2Ae) = 
or 

A = — ^7 and c = 



< 



(31) 



The Ricci scalar of this geometry diverges at r = r s \ := 0, while at r = r S 2 := — e _1 the differential 
scalar V^AV^A diverges. For negative e, the region r > r s i must be removed from the spacetime, 
since otherwise the metric is imaginary, unless c vanishes. For vanishing c and A = 7, we obtain a 
conformally flat solution which may possess one or two horizons surrounding the singularity at the 
origin. For A = —7/2, the spacetime is not conformally flat and may also describe a black hole. 

Let us note that, since in six dimensions the theory is conformally invariant, any metric conformally 
related to (|30p. will be a solution of the system. 



D ^ 6 : For dimensions other than six, the situation is different. Since the theory defined by (|29p 
is not invariant under local Weyl rescaling, one naively expects the factor N (r) in (|28p to be fixed 
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by the field equations. However, this is not the case, and for arbitrary dimensions, the most general 
solution within the family (|28p . for the theory (|29|) is 



ds 2 D = N(r) 



cfcr^ 

(ar 2 + br + j) dt 2 + + r 2 dY? D _ 2 



(32) 



N(r) being an arbitrary function. 

This can be easily seen as follows: Since in dimensions other than six, the trace of the field 
equations for the theory f|29[) is proportional to the Lagrangian, the invariants W\ ~ W2 evaluated 
on a solution should vanish. For the spacetime under consideration (|28p . it has been shown in [22 1, 
that all the components of the Weyl tensor are proportional to a single function X, such that the 
vanishing of X implies that the metric should be conformally flat. Since the restriction W\ = W2 = 0, 
transforms covariantly under Weyl rescalings, it does not involve the function iV (r), and the mentioned 
restriction reduces to X 3 = 0, which implies / (r) = ar 2 + br + 7. Then one is left with a conformally 
flat space, and since the field equations explicitly contain a Weyl tensor, they are fulfilled for any 
arbitrary function N (r For a smooth conformal factor N (r), the spacetime (j32j) is conformally flat 
and it has been studied within the context of conformal gravity in four dimensions in [8] for 7=1 
and in reference |24j for arbitrary 7. The three dimensional cousin of this metric, in which cE 7 is 
replaced by a single compact direction dcfi and 7 is an integration constant, can be obtained through 
"conformal gluing" of BTZ black holes, and is a solution of three dimensional conformal gravity |25| . 
In [IT] this metric was obtained within the context of BHT new massive gravity \10\ [T7] at the special 
point where the two possible maximally symmetric solutions of the theory coincide. In that case, 7 is 
an arbitrary constant, the parameter b plays the role of a gravitational hair, while a is fixed in terms 
of the coupling constant. 

The metric (|32|) may possess an event and a Cauchy horizon, depending on the zeros of gu- It 
generically possesses a curvature singularity located at r = 0, and depending on the sign of the 
integration constant a, it represents an asymptotically locally (A)dS or flat spacetime for (a > 0) 
a < or a = respectively. The details of the different causal structures are given in [TT] , 



4.2 The five-dimensional combination Af e 

As stated in Table I, in five dimensions, there are two linearly independent invariants whose traced 
field equations are of second order. Now, consider the following linear combination as the Lagrangian 

C = l Wx - l -N & (33) 

evaluated on D = 5. This is the unique cubic invariant for which all the components of field equation, 
for static spherically symmetric spacetimes, are of second order. 

n 

7 Note that the same argument is valid for any theory with a Lagrangian of the form C...C provided D 7^ 2n. In 
four dimensions, for n — 3, this solution was found in [23], where it was mentioned that the corresponding model is the 
simplest one that does not admit Schwazschild horizons . 
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As shown in reference |12j . the most general, non-degenerate spherically symmetric solution, is 
given by 

ds 2 = _ C cr 2/3 + \ dt 2 + dr 2 + ^ „ 

\ J _|_ ry I 



where c is an integration constant and 7 = ±1,0 is the curvature of S3. Let us note that this 
spacetime is not conformally flat (it possesses a nonvanishing Weyl tensor) unless c = 0. For positive 
c and 7 = — 1, the metric (I34|) represents a black hole possessing an event horizon located at located 
at r + = c~ 3 / 2 . In this case the geometry of the horizon is given by H3/T, where T is a freely acting 
discrete subgroup of O (3, 1). The horizon hides the curvature singularity located at r = 0, and the 
asymptotic region (r — > 00) is locally flat. Further interesting features of this solution are discussed 
in [T2]. It is also interesting to note that among the class of theories considered here, this is the only 
one which does not admit an (A)dS solution, in the same way as the pure K combination of BHT new 
massive gravity [10J. 

4.3 The three dimensional case 

As shown in Table I, within the family considered, the only nontrivial theory having third order traced 
field equations in three dimensions can be written as C a b c C abc , where C a b c is the Cotton tensor. In 
this theory, the most general static, spherically symmetric solution is given by 



ds 2 = N (r) 



(or 2 + br-n) dt 2 + — 

ar z + or — fi 



dr 2 

2 7/2 

r dm 



(35) 



where a, b and [i are integration constants and N (r) is an arbitrary function. For smooth N(r), 
as mentioned above, this metric has an event and a Cauchy horizon, depending on the value of the 
parameters. 

It will be interesting to study the thermodynamical properties of the black hole within the context 
of AdS/CFT. 



5 Summary 

In this work, we have classified all the six-derivative Lagrangians of gravity for which the trace of 
the field equations have a reduced order. We have seen that, in dimensions greater or equal to six, 
when the trace of the field equations from a generic Lagrangian is restricted to order two, we obtain 
an arbitrary linear combination of three linearly independent curvature invariants, namely, the six- 
dimensional Euler density and the two independent Weyl invariants. These invariants are no longer 
independent in lower dimensions due to the Schouten identities. However, in five dimensions, there is a 
special invariant A/e, which also gives field equations with second order trace. These invariants can be 
used to construct interesting cubic theories of gravity that can serve as toy models for higher derivative 
theories. We have also provided a conjecture regarding all the possible invariants of arbitrary order 
which gives second order traced field equations in any dimensions. In addition, we have obtained 
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the general spherically symmetric solutions for a sub-class of such theories in arbitrary dimensions. 
Our analysis shows that this is possible due to the reduced order of the trace of the field equations. 
When the order of the trace is restricted to three, we obtain two further invariants £ and in 
arbitrary dimensions. These two invariants are not globally independent in three and six dimensions. 
In six dimensions, they reduce to the third type-B anomaly whereas in three dimensions they 
are equivalent to the square of the cotton tensor ~ C abc C abc . We have further obtained a general 
spherically symmetric solution of this theory in three dimensions. 
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A Trace of the field equation 

Here, we prove a general property of the trace of the field equations for any Lagrangian of the form 
C(g a b, Rabcd, V e ) with a fixed degree of differentiation n. It has been shown in [26] that the Lagrangian 
can always be reexpressed as 

£>[9ab> Rabcd, V ai Rbcde, ' ' ' > V( )Rbcde] (36) 

The field equations obtained by variation of the action with respect to the metric takes the form 

_ T ab = _^±_ + E a Rhode + 2 \/ C V d E acdb + -g ab C (37) 
OQab ' 2 

b d ®£* dC, 

dRbcde ^ dV ai Rbcde ^ <9V( Q1 • • • V ap )Rbcde ' 

where T ab is the energy-momentum tensor of the matter fields. Taking the trace of the field equations 
we obtain, 

- T a a = g ab — + E abcd R abcd + -C + tot. deriv. (39) 

OQab 2 



In fact, the anomalies are called global conformal invariants. It was first conjectured by Deser and Schwimmer |28] . 
that any global conformal invariant can be expressed as a linear combination of the Euler density and the local con- 
formal invariants. Recently, the conjecture has been proved by differential geometric techniques by Alexakis [30] and 
cohomological techniques by Boulanger [25]. For six-derivative invariants, in arbitrary dimensions, in addition to the 
two independent Weyl invariants which are purely algebraic, there is a third local conformal invariant which involves 
derivative of the curvature. Our analysis shows that this invariant does not give field equations with reduced order trace 
in arbitrary dimensions. However, they coincide (equivalent upto a total derivative [2S]) with our E and in dimensions 
three and six. 
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Now, if the Lagrangian is of fixed n, then it can be expressed as a linear combination of terms of the 
form 

b"P [R...] q2 [V.i?...P • • • [V. • • • V.R....] qp+2 , (40) 

v ' 

p times 

such that 

2^2 + 3^3 H h (p + 2)q p+2 = n. 

p times 

Then, under the scaling g ab -)• t~ 1 g ab , R bcde -)■ tR bcde , V ttl • • • V ap R bcde -)■ tV ai ■ ■ ■ V ap R bcde , 
the Lagrangian scales as C -> t~ qi+q2+ - q ^+ 2 C. However, gi can be expressed in terms of other g p 's as 

Qi = ~ [4g 2 + 5g 3 + • • • + (P + 4)g p+2 ] 

= ^[2(^2 + 93 + • • • + q P +2) + (2q 2 +3q 3 + --- + {p + 2)q p+2 )} 

= {q2 + qz + ---+q P +2) + -- (42) 

This implies that the Lagrangian scales as t~2£. Now, one can apply Euler's theorem of homogenous 
functions to write the following relation 



— 9 o „ b "i~ Rbcde Q n ^ ai Rbcde Qv -, + • • • + ^(oi ' " ' ^ a„)Rbcde Qv - v— r , 

2 9^" dR bcde dV ai Rbcde v p ^ 3V (oi ■■■V ap) R bcde 

dC dC dC dC 

= 9ab^ 1" RbcdeT-p. 1" V aiJ K;, cde — 1 hV( 

%afe SiZftede dV ai R bcde PJ 5V( ai • • • V ap )Rbcde 

dC 

= g abir - + E bcde R bcde + tot. derm (43) 

Therefore, the trace of the field equations can be written in the form 

n — D 

T\ = — - — C + tot. deriv. (44) 

Now, suppose that the trace of the field equations, from a Lagrangian of n = 6, is of third 
order. Then it must be some linear combination of the invariants L\, ■ ■ ■ , L\ 2 . According to ()44p . in 
dimensions D ^ n, the Lagrangian must be proportional to this combination up to a total derivative. 
In dimensions D = n, since the trace of the field equations is a total derivative, the only way the 
trace can be of at most third order is when it identically vanishes, which is the case for conformally 
invariant theories. 



dC „ dC _ „ dC _ „ dC 
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B Equations of motion 

In this appendix, we provide the details of the analysis for the classification presented in Table I. The 
equations of motion for each term in the general Lagrangian are listed below [27] : 



G iab = ZR aq ef R b qcd Rcdef + 6V p V q (R qcd R p cd ) " (45) 

eg, = m a jR prd R pgr h - m p v q (R p g \R a \ h - R p hbg Rr h ) ~ (46) 

/m(3) p T>cspq r> d p qcd p h n , p dqc jy h p Y7 Y7 f E> P I^P , p p g/ip 

^Sab - K acbdtt Mpqs ~ K a K cdb K qh + ^ b ^adc K qh ~ V p V q {K ah K b + K bh K a H 

+ RCL h R ab + RP h R ab + 2{9 PqR a CdR bhcd + 9ab RWCG lR ' \ cd ~ 9a red ~ 9b R^ 1 ^ red)) 

- \9abL3 (47) 

G fab = ZRapedR^R + R ab RPqCdR pqcd + ^ q{^ RR a b ~ 9a ' 9 b R ' ^ R rscd + 9 Pq 9ab RTSCd R rscd) 

~ \9abU (48) 

G fab = R ac R b fCdR fd + 2R acbd R C f d 9R fg + V p V q (RabR Pq ~ R/R,, 9 + 9 PQ R acbd RCd 

+ g a b RPCqdR cd ~ 9a PR \bd RCd ~ gb'Kad^) ~ \)abH (49) 

G fab = 3R acb dR ec R d + ^V p V q (g pq R a c R bc + g ab R ep R e q - g p R qc R ac - g a p R qc R bc ) - \g a bU (50) 
G fab = R ab RCdR cd + 2RR cd R acbd + V p V q (ga b g pq R cd Rcd + g pq RRab - g a p g b q R cd R C d + g ab RR pq 

- g p RR q - g p RR q ) - ^g ab L 7 (51) 

G fab = 3R2R ab + 3V p V q (g ab g pci R 2 - g p g q R 2 ) - \g ah Ls (52) 

1 
2 

G \oab = V a R ca V b R cd + 2V c R a a V c R bd - U z R ab - V c V d UR cd g ab + V a V c UR c b + V b V c UR c a 

- 2R acbd nR c + 2R c ( a \3R c b * ) + 2V c [/2/V( 6 i? a ) - R d ( a V c R ^ - R dib V a ^R c ] - -g ab Li (54) 
G ?L = " 4 G£l + 2Ggl - + 4Gg> - g£> + 4^ (55) 

r (3) _ r (3) r (3) 1^(3) ,- R . 
^12^ - ^ab ~ '-'Gab + 4^906 ' l bb i 



G 9al = Vai?V 6 i? - 2URR ab - 2(g ab g cd - g ac g bd )V c V d UR - ^g ab L 9 ( 5:! ) 

( 3 ) — T7 Z? crf V7, E> . I 0X7 P d'rjc p. . rn2 p , yy T7 .r~l E> cc ^ 
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Therefore the trace of the full field equations can be expressed as 

^Gg )a = (3 - D/2)A i L i + V p Vg[(6A 1 + 3A 2 - ^-y^A 3 - AA ll )R pabc R q abc 

+ (-3A 2 - 2A 3 + (D - 2)A 5 + 2{D - 2)A 10 + 24A 11 + (D - 2)A 12 )R ab R apbq 

+ {-2A 3 -A 5 + ^ A 6 - 2A 10 + 16A 11 — (D + l)A 12 )R p a R qa 

+ (AA 4 + A 5 + (D- 2) A 7 — (D - 4) A 10 - 12A 11 - ^^-Z^A 12 )RR pq 

+ {-^A 3 + (D- 1)A 4 - A u )g pq R abcd R abcd 

+ (A 5 + ^A e + {D- 1)A 7 - jA 10 - WA 11 - ^A 12 )g pq R ab R ab 

+ (A 7 + 3(D - l)^ 8 - ,4 9 + ^f^A 10 + 3A 11 + ^^-A 12 )g pq R 2 

4 8 

+ (-2(D - l)^ 9 - jA 10 - 2A n - ^-^A 12 )g pq UR] (57) 

Now we impose the trace to be proportional to the Lagrangian. This in turn requires the second term 
on the right hand side to vanish. To realize this, one has to choose the coefficients in such a way 
that, the symmetric tensor quadratic in curvature inside the operator V p V g is proportional to the 
Gauss-Bonnet field equations. This gives us a set of 8 equations in 12 variables and one arbitrary 
parameter u. They are 

6A i + 3A 2 _ ^^A 3 - 4A U = -2u, 

- ^A 3 + (D - l)A* - A 11 = u/2, 

- 3A 2 - 2A 3 + {D- 2)A 5 + 2(D - 2)A W + 24A 11 + (D - 2) A 12 = 4u, 
_ 2A 3 - A 5 + 3{D ~ 2) A 6 - 2A 10 + 16A U - (D + l)A 12 = 4u, 

4A 4 + A n + (D- 2) A 1 — (D - A)A 10 - 12A 11 - ^~^A 12 = -2u, 
+ |a 6 + (D - 1)A 7 - ^A w - WA 11 - l -A 12 = -2u, 

A 7 + 3{D _ i)A S _ A 9 + R^± A 1* + 3A H + R^1 A ^ = u/2, 

4 8 

- 2{D - l)A 9 - ^A 10 - 2A 11 - ^y^A 12 = 0. (58) 

The matrix of linear equations has rank 8, which implies that the general solution can be written in 
terms of 5 arbitrary parameters x,y,z,u and v. In D > 5, the solution is given as 
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A 1 = — [2(D 2 + 5D- W)x + 2(D - 2) 2 y - 6(3D + 2)z + (D 2 - 4)v + 3(D - 2)u] , 

A 2 = -- [8(2D - 3)x + 2(D - 2) 2 y - 8(2D + 3)z + (D 2 - 4)u + 4(L> - l)u] 

A 3 = 2{D - l)x - 2z - u 
A 4 = x 

A" 
A 6 

A 7 
A s 

A 9 



-Ax — (D — 2)y + ^ — --v-2u 
v ,y D-2 2 

1 T 24 

- 8x — 2y — — z + v 

3 y D - 2 



24(£> - l) 2 
1 



A 10 

A 11 
A 12 



4(D - 1)(D - 2) 
4 



[8(D - l)y + 16Dz - (D 2 - D + 2)v - 4(L> - 1) 
[8z + (£> - 2)v] 



D-2 



z — v 



z 

V. 



(59) 



In dimensions D > 5, one can apply the following transformation 

" ~* ~ (0_ 2 )3(0-l) I 4(D " 2)3(D " 1)a + 8(2D " 3)6 + (° 2 + D - 4 H • 

z — > d, 

u^6(D- 5)a, 
w ->■ 24e. 



(60) 



such that 
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where 

£ = — -(3D + 2)Li + -(2D + 3)L 2 - 2L 3 + „ ~ ^ Ls — L 6 ; 2L> xo L 8 

2 V 7 3 V ; D-2 D-2 b 3(D-1) 2 

+ ( J D-2) 2 (D-1) L9 " ^2 Ll ° + Ll1 

1 8D 4L> 4(L>-4) 2D D(D - 3) 

2 v 1 3 £>-2 8 D-2 3(D - l) 2 (D - 2)(D - 1) 

H — — ^~^ 10 + derivative (61) 

G =2(D 2 - 4)Li - 4(D 2 - 4)L 2 - 12DL 5 + 8L 6 + ^_^±^L 8 + — -Lg - 24L 10 + 24L 12 

[D — iy d — i 

/ 9 x / 9 x , s D 2 -D + 2 6D 

=2 [p - 4) L X - 4 (Z? - 4) L 2 - 12(D - 2)L 5 - 16L 6 + — _ L 8 + ^— j-£g 

— 24Lio + total derivative. (62) 

Note that the determinant of the transformation (|60p is . Now in D < 5, one needs to solve 

the system of equations (|5H|) for each value of D separately. The Lagrangians obtained are tabulated 
in Table 1. 
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